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We study the spatially synchronized and temporally periodic solutions of a 1-d lattice of coupled 
sine circle maps. We carry out an analytic stability analysis of this spatially synchronized and 
temporally periodic case and obtain the stability matrix in a neat block diagonal form. We find 
spatially synchronized behaviour over a substantial range of parameter space. We have also extended 
the analysis to higher spatial periods with similar results. Numerical simulations for various temporal 
periods of the synchronized solution, reveal that the entire structure of the Arnold tongues and the 
devil's staircase seen in the case of the single circle map can also be observed for the synchronized 
coupled sine circle map lattice. Our formalism should be useful in the study of spatially periodic 
behaviour in other coupled map lattices. 
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The study of complex dynamical behaviour in extended systems is currently of interest in a wide variety of contexts. 
The modelling and characterization of spatiotemporal behaviour in such systems can provide insights into the 
complex behaviour found in diverse systems like oscillator arrays |^-|^, coupled Josephson junction arrays fl, reaction 
diffusion systems charge density waves [|l2|, biological systems [13|,|4j and turbulence in fluids [[l^ . The 

spatially extended nature of the system permits the appearance of complex spatio-temporal behaviour such as spatially 
periodic, quasi-periodic or chaotic behaviour with the concurrent appearance of temporally periodic, quasi-periodic or 
chaotic behaviour. An understanding of the rich variety of modes which can be excited due to the interplay between 



spatial and temporal behaviour may provide a clue to phenomena like pattern formation in natural systems 17 
and turbulence. 

Coupled Map Lattice (CML) models have recently attracted much attention in the study of spatio-temporal chaos 
and pattern formation as models of spatially extended systems ||^,^. A CML is a dynamical system with discrete 
time, discrete space and continuous states. It consists of dynamical elements on a lattice which are suitably coupled. 
Such systems have been used succesfully to model real life phenomena like spiral waves p6[ | and spatio-temporal 
intermittency jl^ . Though CML models are idealized systems, they are sufflciently complex to be capable of capturing 
the essential features of the dynamics of the system, and at the same time have the advantage of being mathematically 
tractable and computationally efficient. 

Due to the large number of degrees of freedom in such spatially extended systems, a variety of spatio-temporal 
phenomena, like synchronization, intermittency, and spatio-temporal chaos are observed. One of the most important 
and interesting modes which can arise in such systems is the mode corresponding to synchronised behaviour, i.e. 
behaviour in which evolution at each spatial location is identical with that at every other spatial location at any 
arbitrary instant of time. (SeeFig.la) Such synchronized behaviour has been observed in a number of spatially 
extended systems Q, like coupled oscillator arrays [p|,ppC|], coupled pendulums electronic oscillator circuits 
p2| and in pattern formation |2^. As coupled sine circle map systems |^ can constitute models which capture 
many of the essential features of the behaviour of such systems, we address the problem of synchronized or spatially 
homogeneous solutions in a 1-dimensional array of coupled sine circle maps. The temporal evolution of these states 
may be periodic, quasi-periodic or chaotic in nature. The present paper concentrates on synchronized states that are 
temporally periodic in behaviour. 

The single circle map ,is represented by 

et+i^f{et)^et + n--^sM27ret) modi (i) 

where 9t is an angle, at time t which lies between and 1, K is the nonlinearity parameter and Q is the period of the 
system for K — 0. This is one of the simplest representations of physical phenomena involving periodic motion. This 
simple dynamical system which exhibits a tendency to mode lock, as the parameter K is increased, is particularly 
suitable for the description of resonances between periodic motion . An extensive study by Bak et al |^ , shows 
the presence of the Arnold tongues in the H. — K space and complete mode-locking, namely the Devil's staircase 
at K = 1. Since the single sine circle map exhibits this tendency to mode-lock, it is interesting to study whether 
an array of such sine circle maps on a lattice suitably coupled, would demonstrate such a mode-locking spatially. 
Synchronization is the simplest example of such a phenomenon and we focus our attention on the synchronized states 
which are temporally periodic in an array of such coupled circle maps. 

The specific model under study, is a 1-dimensional coupled map lattice of sine-circle maps with nearest neigh- 
bour diffusive symmetric normalized coupling (also called future coupled laplacian coupling) and periodic boundary 
conditions, and is given by 

9t+i{^) = (1 - ^)f{Ot{i)) + y{Ot{i + 1)) + - 1)) modi 

K 



{l-e)[9,{^) + n-—sm{2net{^))) 
K 

W) 



i\et{t + l) + n--^ sin(27r^?t(z + 1)) 



+ et{i-l) + n- -^sm{2T:et{i-l))\ modi (2) 
(27r) J 

where 9t{i) is the angular variable associated with the zth site, at time t. 

The parameters J7 and K are taken to be uniform at each site and defined as in Eq.2 for the single circle map, 
and e which lies between and 1 is the strength of the coupling parameter. 

We investigate in detail the phenomena of spatial synchronization in a system of such coupled sine circle maps. 
We identify the spatially homogeneous and temporally periodic, quasi-periodic and chaotic modes of the system. 
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We carry out a linear stability analysis to analyse the stability properties of the spatially synchronized, temporally 
periodic solutions. The independent variables of the problem are identified and the analysis is cast in terms of these 
independent variables. This leads to a neat and simple block diagonal form for the stability matrix. This permits us 
to identify the regions in parameter space where synchronized solutions of different temporal periods are stable. The 
limits of the mode-locked interval for the temporal period one case can be explicitly evaluated due to the fact that 
the stability matrix has a block diagonal, block circulant structure, for this case. This analysis is also extended to 
higher spatial periods. Our method is quite general and should prove to be useful in the study of spatially periodic 
behaviour in other coupled map lattice models as well. 

The paper is divided into five sections. Section I discusses the coupled shift map, which is the linear version of the 
coupled circle map, namely if = in Eq .2, and is thus the simplest case. We identify the independent variables for 
this system and recast the equation of evolution in terms of these variables. We carry out a systematic linear stability 
analysis, and obtain a neat form for the stability matrix. In section II, we extend the formalism to the system of 
coupled sine circle maps, taking into account the nonlinear sine term as defined in Eq. 2 and cast the problem in 
terms of the new independent variables of this system. The linear stability analysis in terms of these variables, gives a 
neat block diagonal form for the stability matrix. This is followed by an explicit calculation of the eigenvalues for the 
synchronized solution and the evaluation of the limits of the mode-locked interval for the fixed point case. Section III 
presents the generalization to higher spatial periods. In section IV, the numerical simulations for the higher temporal 
periods, of the synchronized solution, are discussed. We extend the algorithm developed by Bak et al |^^, for the 
single circle map, to the synchronized solution of a system of coupled sine circle maps. Section V concludes with results 
and a short discussion. 



I. ANALYSIS FOR THE SHIFT MAP LATTICE 



We begin by carrying out the analysis for the simplest case of Eq. 2, namely for if = 0, which is the case of coupled 
shift maps. 

The single shift map i.e., the linear case of the single circle map, namely if = in Eq.l is given by 

0t+i =Ot + n mod 1 (3) 

This has periodic orbits for rational values of i. e. Q = ^ where P and Q are any integers, and quasi-periodic 
orbits for irrational values of fi. 

Eq. 2, with if = 0, for the coupled shift map, for a lattice of N sites, retaining periodic boundary conditions, 
reduces to 

et+iii)^il-e)(^9tii) + n^+^{etii + l) + 9tii-l))+2n modi (4) 

For the synchronized solution, since the value of the variable at all sites is the same, we note that differences of 
nearest neighbour variable site values is equal to zero for all neighbours. Writing an equation for such a difference for 
the first pair of sites of an array of TV sites, for the coupled shift map, we obtain, 

0i+i(l)-0i+i(2) = (l-e)((?i(l)-0i(2)) 

+ i{9t{2)-9t{3) + 9ta)~9t{N)) (5) 

It is clear that Eq.5 can be expressed completely in terms of the differences defined by 

atii) ^ Otii) ^ 9t{i + I) (6) 
The evolution equation of the differences is given by, 

at+i(i) = (l-e)(at(i)) + |(at(j + l) + modi (7) 

Thus, for the synchronized solution, of the coupled shift map, it is sufficient to define just N differences, to obtain 
the N equations of evolution in terms of independent variables. 

Substituting at(i) = 0, in Eq.7, gives at-^^l{^) = 0, and hence shows that at(i) = 0, is a spatial fixed point solution. 

Expanding up to the linear term about this solution leads to a stability matrix J^^^^* ^ given by 
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(1-6)/ 



(8) 



This is an N X N matrix, which is also circulant and whose eigen values maybe explicitly obtained analytically. 
The eigenvalues of J^^'^f* are given by p6|,p7[ 



A. = (l-e) + -K+c.,^-i) 



(9) 



where LUr is the iVth root of unity given by 



exp(^^) (10) 
On simplifying, this can be written as 



A,^(l-e) + 6Cos( ^^^'^ (11) 

The largest eigenvalue of Eq. 11 namely A''"'^'^^* < 1 defines the stability condition for temporally stable and 
spatially synchronised periodic orbits of the coupled shift map. The largest eigen value is +1, indicating that the 
coupled shift map is marginally stable. We find that as long as the differences between nearest neighbour sites tend 
to zero, we obtain synchronized solutions, irrespective of the corresponding temporal period. Since, the largest eigen 
value does not depend on e, we conclude that, for the synchronized solution, there is no dependence of the coupling 
parameter on the temporal behaviour. Thus for the synchronized solution, the coupled shift map continues to have 
periodic orbits for rational values of f2, and quasi-periodic orbits for 17 irrational, similar to the single shift map and 
there is no enlargement in the phase space due to the coupling. 



II. LINEAR STABILITY ANALYSIS FOR THE COUPLED CIRCLE MAP LATTICE 



We now carry out a similar linear stability analysis with the non-linear terms i.e. for the coupled circle maps as 
defined in Eq.2 with K ^ 0. We consider a 1-d closed chain of N lattice sites, with the sine circle map at each site, 
coupled to its nearest neighbours via diffusive normalized symmetric coupling and periodic boundary conditions, as 
in Eq.2, such that the right hand neighbour of the Nth. point is the first lattice point. 

As discussed in the case of the coupled shift maps, for a synchronised solution, it makes sense to consider the 
evolution of the differences. At any arbitrary time t we write an equation similar to Eq.5, now for the coupled sine 
circle map. Consider, as an example, the first pair of lattice sites, of a lattice of N sites, with evolution at each site 
defined as in Eq.2. 

0t+i(i) - et+i{2) = (1 - e)[(et{i) - et{2)) 

- ^ sin (^(0t(l) - 0t(2))) cos (^((?t(l) + 0t(2))) } 

+ '-{(6,{i)-6m) + {et{2)^et{?,)) 

- ^ sin (^(0t(l) - 0t{N))) cos (^(0t(l) -f 0t{N))) 

~ ^ sin (^(0t(2) - 0t(3))) cos (^((?t(2) + 0t(3))) } mod 1 

(12) 

Eq.l2 shows, that for the coupled sine circle map, the evolution equation for differences between the variable values 
of the nearest neighbours involves not just the differences between pairs of neighbours but also the sum of the variable 
values of the nearest neighbours, unlike the coupled shift map where the evolution equation could be defined in terms 
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of the differences alone. However, it is interesting to note that for a synchronised solution, the difference of nearest 
neighbour lattice sites is zero, and the sum of nearest neighbour lattice sites is a constant, for a fixed fl and K. 
Therefore, we write a second equation for the sums of the variables at nearest neighbour sites, and obtain. 





9t{l) + 0t{2)) 




- - sin (w( 


^et(l) + 0t(2))) cos (n(et{l)-' 


Ot{2)))} 








- - sin (7r( 


'0t(l) + 0t(7V)))cos (71(0^(1)- 


OtiN))) 


- :^sin (n{ 


^0t(2) + ^t(3)))cos (n[0t{2)-, 


614(3))) 1 + 2n modi 



(13) 

which again involves both the sums and the differences of nearest neighbours. It is clear that for the synchronized 
solution of a lattice of coupled sine circle maps, we now need to define equations for the evolution of differences 
of nearest neighbour sites and N equations for the evolution of sums of nearest neighbour sites as compared to N 
differences alone for the coupled shift map. This is due to the fact that the shift map is a linear version of the 
circle map and in the case of coupled circle maps the non-linear sine term requires the identification of the second 
independent variable. At any time t, we define the differences and sums, of the nearest neighbours, as follows, 

atii) = etii)-et{i + i) (14) 
bt{i) = et{i) + et{i + i) (15) 



V i;l,...N 



Using Eq.2, we find that that at{i) and bt{i) evolve via the following Eqs. 

K 



at+i{i) = (1 - e)(at{i) - — sin(7rat(i)) cos{Trbt{i)i)j 

+ ^(at(i + 1) - ^ sin(7ra,(i + 1)) cos(7rfet(i + 1))) 



TT 



+ - {at{i — 1) sm{'irat{i — 1)) cos(7r6t(i — 1)) mod 1 

^ TT 



(16) 



and 



bt+i{i) = (1 - e)(bt{i) - — sm{Trbt{i)) cos(7rat(i))) 

+ Ubt{i + l)- — sm{nbt{i + 1)) cos(7rat(« + 1))) 

+ ^ {bt(i — 1) sin(7r6t(i — 1)) cos(7rat (i — 1))) + 2f2 mod 1 

2 V TT / 

(17) 

For a synchronized solution, as mentioned above, at any time t, 

at{i) = (18) 
bt{i) = constant = s (19) 

Inserting these conditions, Eq.l6 reduces to 

at+i{i)=0 (20) 

implying that a((i) = is a spatial fixed point solution. 
Similarly, Eq.l7 reduces to 
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K K 
bt+\{i) = (1 — e)(s sin(7rs)) + e(s sin(7rs)) + 20. mod 1 



(21) 



which is again a constant (not necessarily the same constant) for a fixed f2 and K . 

Eq.21 is not restricted to temporal evolution of any particular kind, implying that this same equation may be used 
to study temporally periodic, quasiperiodic and chaotic behaviour. As an example, 6t+i(*) = the same constants will 
indicate a temporal fixed point solution and so on. 

Thus &((?') = and hi{i) = constant arc solutions of the equations of evolution. We now perform a Taylor expansion 
up to first order about these solutions to obtain the linear stability matrix Jj which is of order 2N x 2N and of the 
form 



Jt = 



B't 



where 



and 



fesAt{l) e„A(2) 

e„A,(l) £,,^,(2) enAt(i) 

e„At(2) esAtii) 

\enAt{l) 



/e,St(l) e„Bt(2) 

e„Bt(l) e,Bt{2) e„Bt(3) 
e„Bt(2) e,S((3) 



enAt{N)\ 





enAt{N-l) e,At{N)J 



enBt{N)\ 





\enBt{l) 



enBt{N-l) esBt{N)/ 



Here 



Each 



( K \ 
At{i)= (^1 - — cos(7rat(i))cos(7r6t(i))j 



and 



Bt{i) = (l - ^sin(7rat(i))sin(7r6t(i))^ 
Imposing the condition at(i) = 0, bt{i) = s the stability matrix Jt, Eq.22 reduces to 

Jt 



Mt 
Mt 



where each Mt is an N x N circulant matrix given by. 



Mt 



/esAtil) e„ii(2) 
€nAt{l) esA,{2) 
enAt{2) 



En At (TV) \ 






(22) 



(23) 



(24) 



(25) 



(26) 



(27) 



(28) 



(29) 



\e„At(l) ••• esA,{N)/ 
For a synchronised solution all the lattice sites have the same value and hence each bt{i) = s Vi, 1, . . . , A^. Thus each 

K 



Atii) = (l-— cos(7rs)) 



(30) 
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We thus have a block diagonal matrix with circulant blocks. In the stability analysis for other coupled map lattices 
carried out so far, it was necessary to find a similarity transformation to obtain the stability matrix in block diagonal 
form |2^. In this case, we no longer have to find such a transformation, but directly obtain the stability matrix in a 
block diagonal form. 

To study the stability of the homogeneous solution, it is sufficient to obtain the eigenvalues of one of the blocks Mt- 
The eigenvalues of Mt are given by Eq.9 ||2^,^, but are now of the form 

\r = (l-e)(l-Xcos(7rs)) 
+ cjr(^|(l - i^cos(7rs))^ 

where ujr is defined as in Eq.lO 

On simplifying, Eq.31 can be written as 

\r = {l- e)(l - is:cos(7rs)) +e(l - if cos(7rs)) cos( — ^— — 

r = l,2,...,7V (32) 

Now, \\r\ < 1, for all r, ensures the stability of the synchronized solution. The largest eigenvalue \i°-'^9est gj-ossing 
1, is the condition for marginal stability [2^ . 

As mentioned earlier, we have not assumed any temporal periodicity while defining at{i) and bt{i). They have been 
defined at any arbitrary time t. Thus, as long as the condition for synchronization i.e. at{i) = and ht{i) = constant 
is satisfied, this formalism may be used to identify and model any kind of synchronised temporal behavior. 

For the synchronized solution, we can analytically obtain those regions in the ^l — K — t space, for which the spatially 
synchronized (i.e. spatial period 1), and additionally temporal period 1 solution, is stable. 



(31) 



A. The Fixed Point case 



For the temporal fixed point, i.e. at+i{i) = at{i) — and bt+i{i) — bt{i) = s, is also a synchronised solution and 
hence also the spatial fixed point of the system. Using Eq.32, the largest eigenvalue of the stability matrix Jt is given 
by 

A''"'f""* = (l-ifcos(^s)) (33) 
Now, using the fixed point condition for Eq.(21), 

K 

s = s sin(7rs) + 2f7 (34) 

TT 



This gives 



1 ,27rfi, , , 

— arcsm(— — ) (35) 

TT K 



Inserting this value of s and using the condition for stability, Ail, in Eq. 33, we obtain the width of the stability 
interval, 

-("(?))= (if w 

This is the same width, as that obtained by Bak et al for the single circle map. This width is the width of 
that region of f2, for which we obtain stable fixed point solutions. Thus, we find that for the synchronized solution, 
the width of the j interval for the coupled circle map case is identical to that obtained for the single circle map . 
However, in this case, it is the width of the of the j interval for each site, and since the value of at each site is the 
same, for the entire lattice. Thus, for the region of 17, defined by Eq.36, the entire lattice, has synchronized stable 
temporal fixed point solutions. 

For higher temporal periods, the stability matrix Jt, can be obtained in the following way. For a period Q orbit, 
Jt is simply given by. 
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where all the notation is the same as defined earlier. 

Imposing, the condition for a synchronized solution, namely, Eqs.18 and 19, Jt reduces to the following form 

t=i ^ 

The widths of the the higher temporal periods, for the synchronized case cannot be obtained analytically,as has 
previously been observed in the case of the single circle map . To calculate the widths of the higher temporal 
periods, for the synchronized solution we carry out numerical simulations which have been discussed in section IV. A 
detailed discussion is presented in p^ ]. 

Since the formalism set up here is sufficiently general, it may also be used to model synchronisation in other coupled 
map systems. 

The next section extends the formalism for solutions corresponding to arbitrary spatial periods. 



III. ANALYSIS FOR THE KTH SPATIAL PERIOD 



We find that the formalism discussed above may be efficiently used to study the stability of higher spatial periods 
as well. We extend it to higher spatial periods with the stability matrix retaining its block diagonal form. By higher 
spatial periods, say spatial period 4, we mean every ith and (z + 4)th lattice site will have the same value. (See Fig. lb) 

We now extend this formalism, to an arbitrary spatial period k. At any time t, for a spatial period k solution, the 
value of the variable at the ith site is the same as the value at the {i + k)th site. Thus, the difference, between the 
variable values of the ith and the (i + fc)th site, approaches zero, while the sum of these variable values, approaches a 
constant. As discussed in section II, now for a spatial period fc, we can show that the differences and sums, are again 
the independent variables of the system and are now defined as, 

a^{i)^etii)-etii + k) (39) 



b';{i)^et{i) + 9tit + k) (40) 



where the superscript k denotes the spatial period at any time t. 
Eqs.l6 and 17, now for a spatial period fc, are modified to. 



and 



at+i{i) = (1 -e)(«t^(*) - — sin(™f(^))cos(^fef(^))) 

e / K 
+ t: + 1) sin(7rat*^(i + 1)) cos(7r&f (i + 1)) 

+ |(at^(i - 1) - :^sin(7ra,^(z - l))cos(^6f(z - 1))) 



&ti(*) = (1 - e){btii) -^8ininbU^))cosinaU^))) 

+ i (fef + 1) - - sin(^6f (* + 1)) cos(^af (* + 1)) 
+ '^(b^i^ - 1) - - sin(^6f (* - 1)) cos(^a,^(j - 1)) 



jd 1 



2Q mod 1 



(41) 



(42) 



As shown in section II, here too it can be easily shown that, Vi, a^{i) = and 6^(i) = Sk, where si,S2, ■ ■ ■ are all 
constants, and are solutions of the Eq. 41 and 42 for a fixed fl and K. 



8 



For a spatial period k solution, we expand about, a^{i) = and 6^ (i) = constant and obtain the stability matrix 



jk 



(43) 



where 



e„^t^(l) e«4(2) e„4(3) 








Ve„A^(l) 



(44) 



e„A,fc(iV-l) e«At^(iV)/ 



and 



/e.B.'^ll) e„i?,'=(2) 
e„B,'=(l) e.Bf(2) e„Bt'=(3) 



B 



Ik 








\enB!^{l) 
where Cg, e„ are defined by Eq.25 Here, each 

Atii) = (l 



(45) 



enBj^{N-l) esB^{N)J 



K \ 
— cos(7raj (i)) cos(7r6j (i)) ) 



and 



{i) = (l - :|sin(7r4(i))sin(^&,*(z))) 



(46) 



(47) 



Imposing the conditions a^{i) = and 6j («) = Sfe, the stability matrix given by Eq.43 again reduces to a block 



diagonal form 



where each is of the form. 



jk _ 
•Jt — 









(48) 



/esA^il) e„A^{2) 
e„A^(l) e„A,^(2) 
enA^(2) 







\enA^{l) 





e,A,^(iV)y 



(49) 



and each 



(50) 



where, each m goes from 1, . . . , A; See Fig. 1(b) for an illustration of spatial period 4. 

Thus even for higher spatial periods, the stability matrix is in a block diagonal form, Eq.48. Finding the eigen 
values of the matrix M^, and the largest eigen value crossing + 1, gives the region of stability of the period k solution. 
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IV. NUMERICAL SIMULATIONS 



We obtained the width of A|^ri(Y)j , the temporal period 1 case, for the synchronized solution analytically in section 

II. The higher temporal periods, Af2(P/Q), for the synchronized solution, are obtained numerically. 

The widths, Ar2(P/Q), as also stated in section II, are those regions oi il, a.t K = 1, for which, the entire lattice, 
will have temporally stable and spatially synchronized periodic solutions of time period Q. 

We extend the algorithm developed by Bak et al in the following manner. 

For a lattice of N sites we define the following vector notation, 

f{e, n) ^ {h{e, q), h{e, d), . . . , f^ie, n)} (51) 

where i denotes the lattice index and each fi{9, ft) is defined by Eq.2 Each 9 is now a vector of the form 

(6)^ {et{l),et{2),...,9t{N)} (52) 
and the parameter, fi, also a vector, is represented as 

(n) ^ n{2), n{N)} (53) 

ri, in principle may have different values at each site, but in this case has the same value at each site. 
For a 1-d array of coupled sine circle maps, i.e. a multi-dimensional case, the stability criterion is obtained by 
examining the eigen values of the following N x N matrix 



t=i 



9fi 


dfi 


9/1 


det(i) 


oet{2) 


det(N) 


a/2 


dh 


9/2 




det{2) 


det(N) 


dfN 


dfN 


dfN 




det(2) 


det(N) 



(54) 



Let {Xi} be the set of eigenvalues of the matrix St- For a period Q orbit to be stable the eigenvalues of the matrix 
St < 1- The largest eigenvalue crossing 1 defines the marginal stability condition. 

Thus for the higher order temporal periods, we seek the solution to the set of Eqs. 55 and 56, 



Pi9,n)^i9) + p 



the condition for closure, and 



ylarge 



1 



(55) 
(56) 



the condition for marginal stability 

For a general P/Q step, and a lattice of N sites, we start with homogeneous initial conditions and perform the 
iterative technique based on the Newton - Raphson method. 

We define 



gi(0*,f1) = {<?i(l),5i(2) 



where 



and so on. 

We further define. 



_dfl_ 

d9t{l)' d9t{2)'---' d9tiN) 

giil)^ f?i9,n)-i9il))-P 



} 



(57) 
(58) 



g{9,n) 



{9,n^ 
{9,n'^ 



(59) 
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where now, 



For the synchronized solution, it can be shown, that the largest eigenvalue of the matrix, St, namely A'"''^'^^*, is 
equal to each of the N components of 9*2 in Eq.58 | p8| |. 

Therefore, Eqs.55 and 56, may now be expressed as finding {9*,il*) such that 



g{e*,n*) = g* = 

Taylor expanding g* , about the initial point of iteration, g{0^ , il'^) — go, up to the linear order, we obtain, 

9* -90 + AAf 



where 



(60) 



(61) 



(62) 



M 



It is found that, for 



/ dgijl) 

' detii) 



dgi(N) 
992(1) 



, 892 (jV) 
\ d0t{l) 



9*^0 



991(1) 


991(1) 


9gi(l) 


dOtiN) 


ao(i) 


dn(N) 


dgi'iN) 


9gi'(JV) 


dgi(N) 


det(N) 


an(i) 


dn{N} 


992(1) 


992(1) 


992(1) 


det(N) 


ao(i) 


dn{N) 




9g2'(JV) 


dg2(N) 


det(N) 


an(i) 


dn[N) 



(63) 



(64) 



And so as a first approximation. 



-M- 



90 



-M- 



90 



no 



(65) 



(66) 



All derivatives can be derived recursively p8[ | 

We start with homogeneous initial conditions, for the synchronized case, and find, to initiate the iteration, it is 
convenient to locate the superstable point {9, O) = {9s, ds), for all the lattice sites, where now the eigenvalues of the 
matrix St are = 0. All steps have been found to an accuracy of 10~^. 

Using the alogrithm discussed above, at K = 1 we obtain the complete devil's Staircase, of periodic orbits, now in 
the 17 — ^ — e space (see Fig. 2) with all the special features seen in the single circle map. 

We also studied various features of the coupled sine circle map lattice at lower values of K and different Si's. Our 
simulations reveal synchronized quasi-periodic orbits and periodic orbits in the Q — K — e space. We thus obtain the 
entire structure of the Arnold tongues ( see Fig. 3) now with the third coupling dimension e. It is clear from ( see 
Fig. 3 that the three possible routes to chaos seen in the case of the single circle map i.e. via mode-locking alone , 
via quasi-periodicity and mode-locking and via quasi-periodicity to chaos, will also be seen for the coupled sine circle 
map lattice. 

Thus for the synchronized solution, we obtain all the features of the single map, now with a third dimension, in 
the form of the coupling parameter. We also observe that for the synchronised solution all the features as seen in the 
single circle map, are lifted exactly to this extra dimension in the parameter space, namely the coupling parameter e, 
for the entire range, < e < 1. 
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V. CONCLUSIONS 



We have set up a formalism to analyze the phenomenon of synchronization in a lattice of coupled sine circle maps. 
We identify the independent variables, the differences and the sums of variable values at neighbouring sites, at{i) 
and bt{i) for the problem and show that, at{i) — and bt{i) = constant, the synchronised solution, are solutions, of 
the system under consideration. Casting the evolution equations in terms of these independent variables leads to a 
stability matrix which is in a neat block diagonal form. In addition, these blocks have a circulant structure, so we 
have a block diagonal matrix with circulant blocks. The temporal fixed point case, for the synchronised solution, i.e. 
the spatial period one and temporal period one solution can be solved explicitly to obtain the corresponding width of 
the stability interval. We have also demonstrated that this analysis can be easily extended to higher spatial periods. 
For an arbitrary spatial period k, we identify the independent variables, now a^{i) and bf{i) and again show that 
flj (z) — and b^{i) = kconstants, are solutions of the system. A linear stability analysis for this case shows that 
the stability matrix retains its block diagonal form. This analysis may also be easily extended to spatially periodic 
behaviour in other coupled map systems. The analysis is one of the simplest presented so far with the additional 
advantage of being transparent and neat. 

We calculate the widths of the higher temporal periods for the synchronized case, by extending the algorithm 



developed by Baket al |25|. Our numerical simulations reveal that the for the synchronised case, the entire structure 
of the Arnold Tongues and the devil's Staircase as observed for the single circle map is lifted to the chain of coupled 
sine circle maps. 

A study of the parameter basins of attraction, i.e. those regions in the 17 — e, where a particular class of initial 
conditions converges to the synchronized solution, which is the attractor, have also been carried out. These results 
which have been presented elsewhere p8|] showed that these regions were found to be fairly large and showed an 
interesting symmetry about fl — 0.5. Random initial conditions also revealed a finite basin of attraction, which was 
also symmetric about fl — 0.5. The numerical simulations for the random initial conditions showed that there always 
exists a finite region of the Q — e parameter space, at K — 1, for which synchronized solutions will be obtained. The 
temporal period for these synchronized regions was found to be 1. pst . 

Investigation of higher spatial periods reveals extremely interesting structures, results of which will be presented 
elsewhere. Our analysis was carried out for the homogenous case which assumed that the parameter values at 
each lattice site were the same. However the formalism set up is completely general and can be easily extended to 
inhomogeneous cases. This should be useful for the study of lattices where some lattice sites function as pinning sites. 

We hope this work will find useful applications in the study of spatially periodic behaviour in systems like Josephson 
junction arrays and coupled oscillator arrays and in studies of pattern formation, in reaction-diffusion systems and 
charge-density waves. 
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FIGURE CAPTIONS 



FIG.l (a) The synchronised or spatial period 1 solution 
(b) The spatial period 4 solution 

Each symbol in the figure represents a distinct value of the variable of a lattice of 8 sites at any time t. 



FIG. 2 The Devil Staircase for coupled sine circle maps. 
The symmetry about Q = 0.5 is clearly seen. 
The widths are calculated for Af2(^) for 

_P_0111213]_2 
Q ~ 1' 1'2' 3' 3' 4' 4' 5' 5" 

This plot is for 8 lattice sites, but for a synchronized solution, 
identical results are seen for any N. 



FIG. 3 Plot of the Arnold tongues, for the coupled sine circle map lattice. 
The tongues have been plotted for Q = £ starting from 

p _ 11 1 2 
Q 1' 1'2' 3' 3 

This plot is for 8 lattice sites, but identical results are seen for any N, 
for the synchronized solution. 
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